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The supernova remnant Cassiopeia A contains the youngest known neutron star which is also the 
first one for which real time cooling has ever been observed. In order to explain the rapid cooling 
of this neutron star, we first present the fundamental properties of neutron stars that control their 
thermal evolution with emphasis on the neutrino emission processes and neutron/proton superflu- 
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cooling of the neutron star in Cassiopeia A is interpreted as being due to the recent onset of neu- 
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superconductivity determines the observed rapidity of this neutron star's cooling is highlighted. 
This is the first direct evidence that superfluidity and superconductivity occur at supranuclear 
densities within neutron stars. 
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1. Introduction: "neutron stars" do exist 

Neutron stars contain the densest form of cold matter observable in the Universe. Larger 
energy densities are transiently reached in relativistic heavy ion collisions, but the resulting matter 
is extremely "hot". Black holes contain a much denser form of matter, but their interior is not 
observable. Two simple arguments can convince us that such stars, very small and very dense, do 



exist. First, consider the fastest known radio pulsar, "Ter5ad" (PSR J1748-2446ad) []21[], and posit 
that the observed period of its pulses, f = 1.39 ms, is its rotational period. Using causality, that is, 
imposing that the velocity at its equator is smaller than the speed of light c, one then obtains 



2kR „ cP 

2n 



t'equator = = < C =^ /? < — = 65km . (1.1) 



This value of 65 km for the radius R is only a strict upper limit; detailed theoretical models indicate 
radii on the order of 10 km. Secondly, assuming that the star is bound by gravity, we can require 
that the gravitational acceleration agr at the equator is larger than the centrifugal acceleration Ocf 
and obtain 

''sr = ^>^cf = O.^R = ^ or ^>-^ ^ p = __>8xlOi3gcm 3. (1.2) 

Obviously, Newtonian gravity is not accurate in this case, but we can nevertheless conclude that 
the central density of these stars is comparable, or likely larger, than the nuclear density Pnuc — 
2.7 X lO''^ g cm^^. Theoretical models show that densities up to lOpnuc are possibly reachable 



[g6p. In short, a "neutron star" is a gigantic, and compressed, nucleus of the size of a city. 

In what follows, we outline the basic properties of neutron stars relevant for describing their 
thermal evolution with emphasis on the neutrino emission processes and neutron/proton superflu- 
idity/superconductivity. This allows us to present simple analytical models of the cooling of a 
neutron star in order to gain physical insight. The results of these analytical models are comple- 
mented by those of numerical simulations in which full general relativity and the state of the art 
microphysics are employed. Finally, an interpretation of the observed rapid cooling of the neutron 
star in Cassiopeia A as due to the recent onset of neutron superfluidity in it core is proffered. The 
role of proton superconductivity in determining the rapidity of Cas A's cooling is addressed. 

2. Neutron Stars: "pure neutron stars" do not exist 

A "pure neutron star", as originally conceived by Baade & Zwicky and Oppenheimer & 



Volkoff pip cannot really exist. Neutrons in a ball should decay into protons through 

n^p + e^ + Ve- (2.1) 

This decay is possible for free neutrons as m„ > nip + 111^. However, given the densities expected 
within the neutron star interior, the relevant quantities are not the masses, but instead the chemical 
potentials /i, (/ denoting the species) of the participants. The matter is degenerate as typical Fermi 
energies are on the order of tens to hundreds of MeV's, whereas the temperature drops below 
a few MeV within seconds after the birth of the neutron star in a core collapse supernova Jl0|]. 
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Starting with a ball of almost degenerate neutrons, the decay of Eq. (2.1) will generate a Fermi 
sea of protons, electrons and anti-neutrinos. The interaction mean free paths of anti-neutrinos (and 
neutrinos) far exceed the size of the star and these can be assumed to leave the star. Thus, the decay 
will be possible until 

H„=IXp + H,, (2.2) 
where we have neglected /ly- Under this condition, the inverse reaction 

p + e^^n + Ve (2.3) 

also becomes energetically possible. Hence, under equilibrium conditions, which a neutron star 
will reach within a few tens of seconds after its birth, j3 reactions such as Eq. ( ^ ) and (2^) will 



adjust the chemical composition of matter to that characteristic of j3 (or chemical) -equilibrium 



determined by Eq. (2.2). A neutron star is not bom as a "ball of neutrons" which may decay 



according to Eq. ([2.1|), but from the collapse of the iron core of a massive star. It is thus born with 



an excess of protons so that the reaction Eq. (23) initially dominates over Eq. ( |2. 1[ ) in order to 



reduce the proton fraction and it is only when the neutrinos escape that the final cold j8 -equilibrium 



configuration, Eq. ( |2.2[ ), is reached. 

Notice that once /ig > ~ 105 MeV, muons will appear, and be stably present with /x^ = jJ-e- 
The condition for the appearance of muons is fulfilled when the density is slightly above Pnuc- 
Thus, in all processes we describe below, there will always be the possibility to replace electrons 
by muons when the density is large enough to allow for their presence. 

Let us consider simple expressions for the chemical potentials ^ : 

iU« = m„ + ^^^^ + V„ and ^p = mp + ^^^ + Vp (2.4) 



m^c'^ + Pf{e)^c^ = Pf{e)c and = ^^ rn^^c^ + pf{}xYc^ , (2.5) 

where PF{i) is the Fermi momentum of species /, and V„ and Vp are the mean-field energies of n 
and p. For the leptons, Vg and are negligibly small, and we have considered that electrons, but 
not necessarily muons, are ultra-relativistic. With a knowledge of V„ and Vp, the two jS -equilibrium 
relations Eq. ( ^!^ and /x^ = ji^ can be solved. With four /i, 's and two equations, a unique solution 
is obtained by imposing charge neutrality 

np=ne + n^ (2.6) 

and fixing the baryon density 

nB = np+np. (2.7) 

With the «,'s and /^,'s known, one can calculate any thermodynamic potential, in particular the 
pressure P and energy density e = pc^. Varying the value of hb gives us the equation of state (EOS): 
P = P{p). Given an EOS, an integration of the Tolman-Oppenheimer-VoLkoff (TOV) equations of 
hydrostatic equilibrium provides us with a well defined model of a neutron star. 



Relativistic expressions for /x„ and /ip also exist, but are omitted here in tlie interest of simplicity. 
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Figure 1: Schematic illustration of the structure of a neutron star; figure taken from [|38[]. 

The potentials V„ and Vp in Eq. ( ^!4| ) turn out to be rapidly growing functions of density, and 
one can anticipate that eventually reactions such as 

+ A + Vg or/and « + ^r~ + Ve (2.8) 

may produce hyperons. Hyperons can appear, and be stable, once the corresponding j8 -equilibrium 
conditions are satisfied, i.e., /i,, = jX/^ or/and n„ + {j.^ = jx^. At the threshold, where Pf{^) = or 
PfiJL^) = 0, one can expect that |Va.| ^ ?ma and jVj;- 1 ^ m-^- and thus }X\ ~ m^^ and jx^^- ~ m-^-. 
Since wza and mj^- are larger than the nucleon mass by only about 200 MeV these hyperons^ are 
good candidates for an "exotic" form of matter in neutron stars. Along similar lines, the lightest 
mesons, pions and/or kaons, may also appear stably, and can form meson condensates. At even 
larger densities, the ground state of matter is likely to be one of deconfined quarks. All these 
possibilities depend crucially on the strong interactions terms, V,, and Vp. As we will not employ 



them in this chapter, we refer, e.g., to [38] for more details and entries to the original hterature. 



^The Z+ is less favored as its j3 -equilibrium condition is /i£+ = /Xp = /i„ — /ig. Heavier baryon are even less favored, 
but cannot a priori be excluded. 
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Figure |l] illustrates our present understanding (or misunderstanding) of the interior of a neutron 
star, with a black question mark "?" in its densest part. The outer part of the star, its crust, is 
described briefly in the following subsection. 

When only nucleons, plus e's and /x's as implied by charge neutrality and constrained by j3- 
equilibrium, are considered, the EOS can be calculated with much more confidence than in the 
presence of "exotic" forms of matter. For illustrative puposes, we will employ the EOS of Akmal, 
Pandharipande & Ravenhall ^ ("APR" hereafter) in presenting our results. 

Although there is no evidence that any observed "neutron star" or pulsar might actually instead 
be a pure quark star, theory allows this possibility. Such a star would be nearly completely com- 
posed of a mixture of up, down and strange quarks, and would differ from a neutron star in that it 
would be self-bound rather than held together by gravity. 

The reader can find a more detailed presentation and entries to the key literature in [p^]. 



The Neutron Star Crust 

In the outer part of the star, where p < p^uc, a homogeneous liquid of nucleons is mechanically 
unstable (spinodal instablity). Stability is, however, restored by the formation of nuclei, or nuclear 
clusters. At the surface, defined as the layer where P = 0, we expect the presence of an atmosphere, 
but there is the possibility of having a solid surface condensed by a sufficiently strong magnetic 
field [^]. A few meters below the surface, ions are totally ionized by the increasing density (the 
radius of the first Bohr orbital will be larger than the inter-nuclear distance when p ~ lO'^ g cm^^). 
Matter then consists of a gas/liquid of nuclei immersed in a quantum liquid of electrons. When 
p PS 10^ g cm"^. He is of the order of 1 MeV and the electrons become relativistic. From here on, 
the Coulomb correction to the electron gas EOS is negligible and electrons form an almost perfect 
Fermi gas. However, the Coulomb correction to the ion EOS is not negligible. From a gaseous state 
at the surface, ions will progressively go through a liquid state and finally crystallize, at densities 
between 10^ up to ~ lO'*' g cm^^ depending on the temperature (within the range of temperatures 
for which the neutron star is thermally detectable). With growing p, and the accompanying growth 
of He, it is energetically favorable to absorb electrons into nuclei and, hence, the nuclear species 
expected to be present have a neutron fraction strongly growing with p. When p~4 — SxlO^'g 
cm^^ (the exact value depending on the assumed chemical composition), one reaches the neutron 
drip point at which the neutron density is so much larger than that of the proton that not all neutrons 
are bound to nuclei. Matter then consists of a crystal of nuclei immersed in a perfect Fermi gas of 
electrons and a quantum liquid of dripped neutrons. This region is usually called the inner crust. In 
most of this inner crust, the dripped neutrons are predicted to form a superfluid (in a spin-singlet, 
zero orbital angular momentum, state ^Sq). All neutron stars we observe are rotating; a superfluid 
cannot undergo rigid body rotation, but it can simulate it by forming an array of vortices (in the 



core of which superfluidity is destroyed). (See, e.g., []48|].) The resulting structure is illustrated in 
the inset A of Fig. |l} 

At not too high densities, nucleons are correlated at short distances by the strong interaction 
and anti-correlated at larger distances by the Coulomb repulsion between the nuclei, the former 
producing spherical nuclei and the latter resulting in the crystallization of the matter. As p ap- 
proaches Pnuc the shape of the "nuclei" can undergo drastic changes: the nuclear attraction and 
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Coulomb repulsion length-scales become comparable and the system is "frustrated". From spheri- 
cal shapes, nuclei are expected to become elongated ("spaghettis"), form 2D structures ("lasagnas"), 
always surrounded by the neutron gas/superfluid which occupies an increasingly small portion of 
the volume. Then the geometry is inverted, with the dripped neutrons confined into 2D, then ID 
("anti-spaghettis") and finally OD ("swiss cheese") bubbles. The homogeneous phase, i.e. the core 
of the star, is reached when p = Pcc — 0.6pnuc ^- This "pasta" regime is illustrated in the inset B of 
Fig. [l] and is thought to resemble a liquid crystal [pO[]. 



3. Neutrino Emission Processes 

The thermal evolution of neutron stars with ages ~ 10^ yrs is driven by neutrino emission. 



Here we describe the dominant processes. The simplest neutrino emitting processes, Eq. (gjj) and 

(3.1) 



Eq. (|2.3[), are 

j n — > p + e^ + V 



p + e — > n + v, 



e 1 



and are generally referred to as the direct Urea ("DU") cycle. 

By the condition of j8 -equilibrium the cycle naturally satisfies energy conservation, but mo- 
mentum conservation is much more delicate. Due to the high degeneracy, all participating particles 
have momenta p{i) equal (within a small T <^Tf correction) to their Fermi momenta As 

1 /3 

PF{i) °^ n- and rip ~ «^ ^ «„, momentum conservation is not a priori guaranteed. It is easy to 
see that, in the absence of muons and hence with rip = rie, the "triangle rule" for momentum con- 
servation requires that Xp > l/9~ 11%, whereas at p ~ Pnuc we have Xp ~ 5%. In the presence of 
muons, which appear just above Pnuo the condition is stronger and one needs Xp larger than about 
15% [27]. The proton fraction Xp grows with density, its growth being directly determined by the 



growth of the nuclear symmetry energy, so that the critical proton fraction for the DU process is 



likely reached at some supra-nuclear density []27[]. For the EOS of APR, the corresponding critical 



neutron star mass for the allowance of the DU process is 1 .91Mq , but other EOS 's predict smaller 
values. 

At densities below the threshold density for the DU process, a variant of this process, the 
modified Urea (MU) process 

^P + n' + e- + v, (3.2) 
— > n + n' + Ve 




can operate as advantage is taken of neighboring particles in the medium []17p. A bystander neutron 
n' can take or give the needed momentum for momentum conservation. The processes in Eq. (3.2) 
show the neutron branch of the MU process and in the proton branch n' is replaced by a proton p'. 
As it involves the participation of five degenerate particles, the MU process is much less efficient 
than the DU process. Unlike the DU processes which require sufficient amount of protons, both 
branches of the MU process operate at any density when neutrons and protons are present. 



Pnuc correspond to the density of symmetric nuclear matter, i.e. with a proton fraction Xp = 50%, and zero pressure, 
whereas in a neutron star crust at p ~ p^uc one has ~ 3 — 5%. 
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Name 



Process 



Modified Urea 
(neutron brancli) 
Modified Urea 
(proton brancli) 

Bremsstralilungs 

Cooper pair 

Direct Urea 
(nucleons) 
Direct Urea 
(A iiyperons) 
Direct Urea 
(Z^ hyperons) 
71^ condensate 
K condensate 



Direct Urea cycle 
(u-d quarks) 



Emissivity 
(erg cm^^ s^^) 



n + n^n + p + e^ + Ve 
n + p + e^ n + n + Ve 
p + n—^p + p + e^ + Ve 
p + p + e^^p + n + Ve 
n + n -^n + n + v + V 
n+p^n+p+v+v 
p + p p + p + V + V 
n + n^ [nn] + V + V 
P + P^[pp] + V + V 
n-^ p + e^ + Ve 
p + e^ ^n + Ve 
p + e^ + Ve 
p + e^ ^ K+Ve 

n + e^ + Ve 

n+ < 71" >-^ n + e^ + Ve 
n+ < K >-^ n + e^ + Ve 



1021 R 



- 10^"^ RT^ 

5x10^^ RT^ 
5x10^''/? To'' 



10^^ RT^ 



IQ^'^RT^ 



10^'^ RT^ 

lO^^RT^ 
10^^ R 



d ^ u + e + Vg 
u + e^ — 7- (i + Vg 



10 



27 



RTl 



Slow 
Slow 

Slow 

Medium 
Fast 
Fast 

Fast 

Fast 
Fast 



Fast 



Direct Urea cycle 


5 — )■ M + e + Vg 


~ io2^/?r| 


Fast 


(u-s quarks) 


u + e^ — )■ 5 + Vg 







Table 1: A sample of neutrino emission processes. is temperature T in units of 10^ K and the /?'s are 
control factors to include the suppressing effects of pairing (see §^. 



In the presence of hyperons, DU and MU processes which are obvious generalizations of 



the nucleon-only process can also occur [42]. When they appear, the A's have a density much 
smaller than that of the neutron and hence a smaller Fermi momentum. Consequently, momentum 
conservation in the DU cycle p + e^ — A + Vg and A — + e" + Vg is easily satisfied, requiring 
an of only ~ 3%. Notice that if the nucleon DU process is kinematically forbidden, the hyperon 
DU process £" — )• ?i + e" + Vg with n + e^ — £" + Vg is also kinematically forbidden, whereas the 
no-nucleon DU process £" — )- A + e" + Vg together with A + e" — ^ £" + Vg may be possible and 
require very low threshold fractions Xf^ and . 

In deconfined quark matter, DU processes such as u + e ^ d + Ve and J — )• m + e + Vg are 



possible IgSp. Reactions in which the d quark is replaced by an s quark can also occur in the case s 
quarks are present. 

In the presence of a meson condensate copious neutrino emission in processes such as 



n + e +Vg 
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and 

— >n + K^ + Ve 




(3.4) 

n + e +Ve 



occur As the meson condensate is a macroscopic object there is no restriction arising from 

momentum conservation in these processes. 

Finally, another class of processes, bremsstrahlung, is possible through neutral currents JI^]. 
Reactions such as 

n + n — >n + n' + Ve + 'Ve (3.5) 

P + P > p + p +Ve + Ve (3.6) 

n + p — >n + p + Ve + Ve (3.7) 

are less efficient, by about 2 orders of magnitude, than the MU processes, but may make some 
contribution in the case that the MU process is suppressed by pairing of neutrons or protons as we 
will see in §|5[ 

In Table | we list these processes with order of magnitude estimates of their neutrino emissiv- 
ities. Most noticeable is the clear distinction between processes involving 5 degenerate fermions 
with a r*^ dependence, which are labeled as "slow", and those with only 3 degenerate fermions with 
a dependence, which are several orders of magnitude more efficient and labeled as "fast". The 
difference in the T dependence is important and simply related to phase space arguments which are 
outlined below. The Cooper pair process [|15, 51] will be described in §|5[ 



The reader can find a detailed description of neutrino emission processes in [ 54 ] and an alter- 
native point of view in [^]. 

Temperature dependence of neutrino emission 

We turn now to briefly describe how the specific temperature dependence of the neutrino pro- 
cesses described above emerges. Consider first the simple case of the neutron beta decay. The 
weak interaction is described by the Hamiltonian Ml = {Gp / V^)Bp^L^ , where Gp is Fermi's 
constant, and L>^ = \j7^Y'^{l - Y5)Yv and = \j7^Y^{Cvl — Ca75)v^« are the lepton and baryon 
currents, respectively. In the non relativistic approximation, one has = cos dc^J.,^^ and B' = 
— cos dcgA ^],<y'^n where dc is the Cabibbo angle and the axial-vector coupling. Fermi's 
Golden rule gives us for the neutron decay rate 

i.e., a sum of {2nY5'^{Pf — Pi) ■ \Mfi\^ over the phase space of all final states / = {py,pe,Pp). 
The integration gives the well known result: Wp = [G| cos^ 0f ( 1 + 3g^ )mlc^]wp , where w^j ~ 1 
takes into account small Coulomb corrections. This gives us the neutron mean life, T„ ~ 15 min., 
or, measuring T„, a measurement of Gp (modulo cos dc and wp). 

The emissivity of the DU process (the Feynman diagram for this process is shown in 
Fig. |2|) can be obtained by the same method as above leading to 
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Direct Urea Modified Urea Bremsstrahlung 



Figure 2: Feynman diagrams for the indicated neutrino emitting processes. 



with an extra factor for the neutrino energy and the phase space sum now includes the initial 
n. The fi terms, fi being the Fermi-Dirac distribution for particle / at temperature T, take into 
account: (1) the probability to have a « in the initial state, /„, and (2) the probabilities to have 
available states for the final e and p, denoted by (1 — /g) and (1 — /p), respectively. We do not 
introduce a Pauli blocking factor (1 — /v) for the anti-neutrino as it is assumed to be able to freely 
leave the star (i.e., fy = 0). When performing the phase space integrals, each degenerate fermion 
gives us a factor T , as particles are restricted to be within a shell of thickness ksT of their respective 
Fermi surfaces. The anti-neutrino phase space gives a factor T^. The factors Ey and the delta 
function 8'^{Pf — Pi) gives a factor T^^ (it cancels one of the T's from the degenerate fermions). 
Altogether, we find that 

ocT'' -T -T -T ■j-{\f -T = T^ , (3.10) 
where the (1)^ factor emphasizes that the squared matrix element |M/;p is T-independent. An 



explicit expression for the neutrino emissivity for the DU process can be found in [ 27 ] . 

Figure Fig. § shows a Feynman diagram for the MU process. There are two more such di- 
agrams in which the weak interaction vertex is attached to one of the two incoming legs. In this 
case, the T-power counting gives 

ocT^ -T -T -T -T -T -T = T^ , (3.11) 

where the |My;p now involves a strong interaction vertex, the wavy line in Fig. ^, but is still T- 
independent. Notice that in the MU case, the internal neutron is off-shell by an amount ~ pLg which 
does not introduce any extra T-dependence as we are working in the case 100 MeV>r. 



Reference [17] contains the expression from which neutrino emissivity from the MU process can 
be calculated. Considering finally the n — n bremsstrahlung process. One diagram is shown in 
Fig. ^ and there are three more diagrams with the weak interaction vertices attached to the other 
three external lines. The T-power counting now gives 

e^'' ocT^ -T^ -T -T -T -T -j-l^ -7 = 7^ (3.12) 
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with two factors for the neutrino pair and a (T^^)^ from the matrix element: the intermediate 
neutron is almost on-shell, with an energy deficit ~ T, and its propagator gives us a T^' dependence 
for Mfi. A working expression for the bremsstrahlung process can be found in [[T 



4. Neutron Star Cooling 

The basic features of the cooling of a neutron star are best illustrated by considering the energy 
balance equation for the star in its Newtonian formulation"* : 



dt ^ dt 



-Ly — Ly + H , 



(4.1) 



where Ei^ is the thermal energy content of the star, T its internal temperature, and Cy its total 
specific heat. The two energy sinks are neutrino luminosity Ly, described in §^, and the surface 
photon luminosity Ly, discussed in § 1^. The source term H would include heating mechanisms as. 



e.g., magnetic field decay, which we will not consider here. 
4.1 Specific Heat 

The dominant contributions to the specific heat Cy come from the core, constituting more than 
90% of the total volume, whose constituents are quantum liquids of leptons, baryons, mesons, and 
possibly deconfined quarks at the highest densities. Hence one has 



Cv = 52^v,i with Cv, 



Cyjdv {i = e, n, n, p, hyperons, quarks, 



(4.2) 



where Cy,, is the specific heat, per unit volume, of component /. For normal (i.e., unpaired) degen- 
erate fermions, the standard Fermi liquid result [^] 



N{Q)—k%T with A^(0) 



(4.3) 



E 
o 



0) 

o 

CM 
O 




1 

Stellar Volume [10'° cm"*] 



Figure 3: Cumulative specific heats of e, jU, p, and n vs. stellar volume in the core of a 1.4 Mq star built 
using the EOS of APR, at temperature T = 10^ K. Nucleons are assumed to be unpaired. No hyperons or 
quarks are present in the EOS. This figure is adapted from [^5[]. 

^Numerical results to be shown later include full general relativistic effects. 
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can be used, where m* is the fermion's effective mass. In Fig. |3|, the various contributions to Cy 
are illustrated. When baryons, and quarks, become paired, as briefly described in §§, their con- 
tribution to Cv is strongly suppressed at temperatures T <^Tc (Tc being the corresponding critical 
temperature). Extensive baryon, and quark, pairing can thus significantly reduce Cy, but by at most 
a factor of order ten as the leptons do not pair. The crustal contribution is in principle dominated by 
neutrons in the inner crust but, as these are certainly extensively paired, practically only the nuclear 
lattice and electrons contribute. 



4.2 Photon emission and the neutron star envelope 

Thermal photons from the neutron star surface are emitted at the photosphere, which is usually 
in an atmosphere but could be a solid surface in the presence of a very strong magnetic field. The 
atmosphere, which is only a few centimeters thick, presents a temperature gradient. Consequently, 
photons of increasing energy escape from deeper and hotter layers. It is customary to define an 
effective temperature, Te, so that the total surface photon luminosity, by analogy with the blackbody 
emission, is written as 

Ly = 4KR^OsBT,\ (4.4) 

where ctsb is the Stefan-Boltzmann constant. Observationally, Ly and are red-shifted and Eq. ( p!4| ) 
is rewritten as 

Ly„ = 47iRlasBT,t (4.5) 

where Lyoo = &^^Ly, Tgoo = e'^T^, and = e^'^R. Here e^'^ = 1 + z, with z being the red-shift, and 
e^'^ is the ^oo coefficient of the Schwarzschild metric, i.e.. 



2GM 

Notice that R^ has the physical interpretation of being the star's radius that an observer would 
measure trigonometrically, if this were possible. 

In a detailed cooling calculation, the time evolution of the whole temperature profile in the star 
is followed. However, the uppermost layers have a thermal time-scale much shorter than the interior 
of the star and are practically always in a steady state. It is hence common to treat these layers 
separately as an envelope. Encompassing a density range from ph at its bottom (typically ph = 10^" 
g cm^^) up to that at the photosphere, and a temperature range from Th to Te, the envelope is about 
one hundred meters thick. Due to the high thermal conductivity of degenerate matter, stars older 
than a few decades have an almost uniform internal temperature except within the envelope which 
acts as a thermal blanket insulating the hot interior from the colder surface. A simple relationship 



between Tj, and Tg can be formulated as in 1 19] 



(rp \ 0.5+a 
Wk) 

with a ^ 1 . The precise Tg — Th relationship depends on the chemical composition of the envelope. 
The presence of light elements, resulting in larger thermal conductivities, implies a larger Tg for the 
same Th compared to the case of a heavy element envelope. Magnetic fields also alter this T^ — Th 



relationship (see, e.g., [34] for more details) 
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4.3 Some simple analytical solutions 

As the essential ingredients entering Eq. ( p~T| ) can all be approximated by power-law functions, 
one can obtain simple and illustrative analytical solutions. Let us write 

Cy = C9-T<), Ly=Ng-T^ and Ly = S9-T^+^", (4.8) 

where Tg = r/(10^ K). As written, Ly considers slow neutrino emission involving five degenerate 
fermions from the modified Urea and the similar bremsstrahlung processes, summarized in Table |l|. 
The photon luminosity Ly is obtained from Eq. ( p!4| ) using the simple expression in Eq. (p77|). 
Typical values are C9 ~ 10^^ erg K^^,N<) ~ 10^" erg s^^ and ^9 ~ 10^^ erg s^' (see Table 3 in [34] 
for more details). In young stars neutrinos dominate the energy losses (in the so-called neutrino 
cooling era) and photons take over after about 10^ yrs (in the photon cooling era). 
Neutrino cooling era: In this case, we can neglect Ly in Eq. ( p~T| ) and find 



WC9 ( 1 



r = IO'^K(tmuA)'^^ (when r< To) (4.9) 



with a MU cooling time-scale Tmu = 19^C9/6A'^9 ~ 1 yr when the star reaches the asymptotic 

solution (T <^Tq,Tq being the initial temperature at time t = to = 0). 

Photon cooling era: In this era, Ly becomes negligible compared to Ly so that we obtain 

'-" + 4«V0-^^>° ^ (whe„T«r,) (410) 

where Ti is T at time t\. Notice here that the slope, — l/4a, of the asymptotic solution is strongly 
affected by any small change in the envelope structure, Eq. ( |4.7[ ). 

4.4 Some numerical solutions 

Numerical simulations of a cooling neutron star use an evolution code in which the energy 
balance and energy transport equations in their fully general relativistic forms are solved, usually 
assuming spherical symmetry and with a numerical radial grid of several hundred zones ^. A set of 
cooling curves that illustrate the difference between cooling driven by the modified Urea and the 
direct Urea processes is presented in Fig. |[ Cooling curves of eight different stars of increasing 
mass are shown, using an equation of state model, from [pl|], which allows the DU process at 
densities above 1.25 x lO'^ g cm^^, i.e., above a critical neutron star mass of 1.35 Mq. Notice 
that the equation of state used is a parametric one and its parameters were specifically adjusted to 
obtain a critical mass of 1.35 Mq which falls within the expected range of isolated neutron star 
masses; other equations of state can result in very different critical masses. The difference arising 
from slow and fijst neutrino processes is clear. 

Four successive cooling stages are marked in the figure. The neutrino cooling era is marked 
stage 3 and the photon cooling era is marked stage 4. As the figure shows T^, and since ~ 
r"^, the slopes of the cooling curves in this figure are rescaled from the values of Eq. ( |4.9| ) and 
( |4.10 ), i.e., they become — 1/12 and f» — l/8a, respectively. The simple analytical model of the 



Such a code, NS Cool, is available at: http : / /www .astrosou.unam.mx/neut rones /NS Cool/ 
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Figure 4: Slow versus fast cooling (see text for details of stages 1 through 4 marked in the figure). Figure 
adapted from p3| ] . 

previous section considered a single temperature in the stellar interior. For a young star, whose age 
is smaller that its thermal relaxation time, the interior temperature has a very complicated radial 
profile. In particular, the core cools much faster than the crust, due to its much stronger neutrino 
emission. During stage 1 in Fig. ^, the surface temperature is controlled by processes occurring 
in the outer layers of the crust and is totally independent of the temperature deeper in the star. 
(One can appreciate that the more massive stars have a larger at that time; this is mostly due to 
the fact that their blanketing envelopes are thinner than in low mass stars.) The rapid drop in 
occurring during stage 2 corresponds to the thermal relaxation of the crust; the star's age becomes 
comparable to the heat transport time-scale from the crust to the core and the crustal heat flows 
into the core. After stage 2, the stellar interior is essentially isothermal with a strong temperature 
gradient present only in the envelope and the simple analytical solutions presented above apply. 



5. Pairing and its Effects 

Pairing, which induces superfluidity in the case of degenerate neutral fermions and supercon- 
ductivity for charged fermions, is expected to occur between neutrons/protons in the interior of 



neutron stars [30]. The Cooper Theorem |12] states that, in a system of degenerate fermions the 
Fermi surface is unstable, at T = 0, due to the formation of Cooper pairs if there is an attractive 
interaction in some spin-angular momentum channel. The essence of the BCS theory [||] is that as 
a result of this instability there is a collective reorganization of the particles at energies around the 
Fermi energy and the appearance of a gap in the quasi-particle spectrum (see Fig. ^ which is the 
binding energy of a Cooper pair. At high enough temperature the gap A(r) vanishes and the system 
is in the normal state. The transition to the superfluid/superconducting state is a second order phase 
transition and the gap A(r) is its order parameter (see Fig. Explicitly, A(r) = when T > T^- 
and, when T drops below T^-, A(r) grows rapidly but continuously, with a discontinuity in its slope 
at r = Tc. This is in sharp contrast with a first order phase transition, in which the transition occurs 
entirely at Tc (see left panel of Fig. ^ and will be of paramount importance for our purpose. In the 
BCS theory, which remains approximately valid for nucleons, the relationship between the gap and 
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Normal Fermi Liquid Superfluid Fermions 




Figure 5: Comparison of quasi-particle spectra, e{k), i.e., the pole in the 2-point Green's function, for a 
normal and a superfluid Fermi liquid. The reorganization of particles at e ^ £/r into Cooper pairs results in 
the development a gap 2A{k) in the spectrum so that no particle can have an energy between Ef — A(fcf ) and 
eF+A{kf). 

Tc is 

^{T = Q)c^l.l5kBTc. (5.1) 

In a normal Fermi system at T = 0, all particles are in states with energies £ < £f . When 
r > 0, states with energies £ ~ £f can be smoothly occupied (left panel of Fig. |5|) resulting in 
a smearing of the particle distribution around ep in a range ~ ksT . It is precisely this smooth 
smearing of energies around ep which produces the linear T dependence of Cy, § |4.1[ , and the or 
dependence of the neutrino emissivities, § ^. 

In a superfluid/superconducting Fermi system at T = 0, all particles are in states with energies 
£ < £f (actually, £ < £f — A). When <T <Tc, states with energy e> Ep (actually, £ > £f + A) 
can be populated. However, in contrast to the smooth filling of levels above Sp in the case of a 



First Order Phase Transition Second Order Phase Transition Smooth State Transition 




Tc T Tc T T 



Figure 6: Temperature evolution of the state of a system parametrized by an "order" parameter, A{T). 

First order phase transition: discontinuous change of A at Tc; latent heat due to the entropy difference 

between the two states. (Examples: solid o liquid; liquid gas below the critical point.) 

Second order phase transition: continuous change of A but with a discontinuity in the slope at Tc; no 

latent heat but a discontinuity in specific heat. (Examples: superfluid o normal fluid; ferromagnetic O 

paramagnetic.) 

Smooth state transition: continuous change of A with no critical temperature; increase of specific heat in 
the regime where A(r) changes rapidly. (Examples: liquid gas above the critical point; atomic gas o 
plasma.) 
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Figure 7: Control functions for MU neutrino emission (left panel) and specific heat (right panel). 

normal Fermi liquid, the presence of the 2A{T) gap in the spectrum implies that the occupation 
probability is strongly suppressed by a Boltzmann factor ~ eKp[—2A{T)/kBT]. Consequently, all 
physical properties/processes depending on thermally excited particles, such as the specific heat 
and the neutrino emission processes described in § ^, are strongly suppressed when T <^Tc. In 
practice, for numerical simulations of neutron star cooling, these suppressions are implemented 
through "control functions" such that 



„Paired 



RcC^ 



Normal 



, Paired 



Normal 



(5.2) 
(5.3) 



There is a large family of such functions, for each process and they moreover depend on how many 
of the participating particles are paired and on the specific kind of pairing; a few examples are 
displayed in Fig.^Q. 

5.1 Theoretical expectations on pairing gaps 

Soon after the development of the BCS theory, Bohr, Mottelson & Pine |j8|] pointed out that 
excitation energies of nuclei show a gap, as shown in the left panel of Fig. |l Even-even nuclei 
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Figure 8: Left panel: lowest excitation levels of nuclei (adapted from []8[]). Right panel: binding energy per 
nucleon for the most beta-stable isobars. 
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Figure 9: Left panel: possible spin-angular momentum combinations for Cooper-pairs. Right panel: phase 
shifts for N-N scattering as a function of the laboratory energy (middle axis) or the neutron Fermi energy 
and density for a neutron star interior (lower axis). Adapted from [p7|]. 



clearly require a finite minimal energy for excitation. This energy was interpreted as being the 
binding energy of the Cooper pair which must break to produce an excitation. In contrast, odd 
nuclei do not show such a gap, and this is due to the fact that they have one nucleon, neutron or 
proton, which is not paired and can be easily excited. The right panel of Fig. |8| shows that pairing 
also manifests itself in the binding energies, even-even nuclei being slightly more bound than odd 
nuclei^. 

As a two-particle bound state, the Cooper pair can appear in many spin-orbital angular mo- 
mentum states (see left panel of Fig. ^). In terrestrial superconducting metals, the Cooper pairs are 
generally in the ^Sq channel, i.e., spin-singlets with L = orbital angular momentum, whereas in 
liquid ^He they are in spin-triplet states. What can we expect in a neutron star ? In the right panel 
of Fig. 1^ we adapt a figure from one of the first works to study neutron pairing in the neutron star 
core [ |47[ ] showing laboratory measured phase-shifts from N-N scattering. A positive phase-shift 
implies an attractive interaction. From this figure, one can expect that nucleons could pair in a 
spin-singlet state, ^Sq, at low densities, whereas a spin-triplet, ^P2, pairing should appear at higher 
densities. We emphasize that this is only a presumption as medium effects can strongly affect 
particle interactions. 

A simple model can illustrate the difficulty in calculating pairing gaps. Consider a dilute Fermi 
gas with a weak, attractive, interaction potential U . The interaction is then entirely described by 
the corresponding scattering length, a, ~' which is negative for an attractive potential. In this case, 

"Notice that, as a result of pairing, the only stable odd-odd nuclei are 2h(1,1), *'Li(3,3), "'B(5,5), and ''♦nc?,?). All 
heavier odd-odd nuclei are beta-unstable and decay into an even-even nucleus. 

^The scattering length a is related to C by a = {m/47th)Ui) with t/k = / rf-'r exp((k- r) t/(r). 
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Figure 10: Three collections of predicted pairing critical temperatures for neutrons in the 'Sq and ^Pi 
channels and protons in the ^So channel. See for references to the original works. 

the model has a single dimensionless parameter, \a\kf, and the dilute gas corresponds to \a\kf <^ I. 
Assuming the pairing interaction is just the bare interaction U (which is, improperly, called the BCS 
approximation), the gap equation at T = can be solved analytically, giving the weak-coupling 
BCS-approximation gap: 



8 fh^k^ 

Aikp) ABCsikp) = -2 ( ^ ) exp 



71 



2\a\kf 



when \a\kf — > . (5.4) 



This result is bad news: the gap depends exponentially on the pairing potential U. The Cooper pairs 
have a size of the order of ~ hvp /A (the coherence length) and thus t,kF ^ exp[+7r/2|a'|^f ] » 1 
in the weak coupling limit. There appear to be many other particles within the pair's coherence 
length. These particles will react, and can screen or un-screen, the interaction. Including this 
medium polarization on the pairing is called beyond BCS, and in the weak coupling limit its effect 
has been calculated analyticaly JTs]], giving 



A{kF) — y AGMBikp) = j^^yj^ABcsikp) ^0.45 ABcsikp) when \a\kF->-0. (5.5) 

So, screening by the medium reduces the gap by more than a factor two, even in an extremely dilute 
system. 

Pairing correlations in nuclei are part of everyday nuclear physics, and a significant amount of 



work has also been devoted to the neutron star environment (see, e.g., [ |13| ] and [ |28| ] for reviews). 
We show in Fig. [l^ three sets of predicted Tc for the neutron star interior. At low density, corre- 
sponding to the neutron star crust in the regime of dripped neutrons, the expectation of a neutron 
^So superfluid is amply confirmed by the models. This regime was also illustrated in the inset A 
of Fig. [T[ At higher densities, corresponding to the neutron star core, the situation is much more 
ambiguous. Due to their low concentrations, protons have small Fermi momenta in the core and 
are expected to form a ^Sq superconductor. There is, however, a significant uncertainty in the 
size of their gap, with predicted values of Tc ranging from ~ 10^ K up to 6 x 10^ K, and a larger 
uncertainty in the range of Fermi momenta in which A{kf) is non-zero, which translates into an 
uncertainty of a factor of more than 3 on the density range covered by the superconductor. In the 
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"pessimistic" case protons would be superconducting only in the outer part of the core, whereas in 
the "optimistic" case the whole core may be superconducting. 

For neutrons in the neutron star core, there is no agreement between the many published mod- 
els about either the maximum value of Tc or on the density range in which pairing is significant. 
Notice that, due to the tensor interaction, pairing is expected to be in the ^P2-^F2 channel. However, 
even the best models of the N-N interaction in vacuum fail to reproduce the measured phase shift in 
the channel [Q]. Due to medium polarization a long-wavelength tensor force appears that is not 
present in the in vacuum interaction and results in a strong suppression of the gap [p3|]. Recently, 
the effect of three-body forces (TBF), absent in the laboratory N-N scattering experiment, has been 
considered. TBF are necessary to reproduce the nuclear saturation density; they are, in the bulk, 
repulsive and their importance grows with increasing density. However, it was found in [57, 58] 
that, at the Fermi surface, they are strongly attractive in the ^P2-^F2 channel and result in very large 
neutron ^P2-^F2 gaps. Other delicate issues are the effect of the proton contaminant and the likely 
development of a tt" condensate^ which also strongly affects the size of the neutron (and proton) 
gap(s). In short, the size and extent in density of the neutron ^P2-^F2 gap in the neutron star core is 
poorly known. 



5.2 The Cooper pair neutrino process 

The formation of the fermonic pair condensate also triggers a new neutrino emission process. 



which has been termed as the "pair breaking and formation", or PBF, process ||44|]. When an ff pair 
(f = n, p, or any fermion undergoing pair condensation) forms, its binding energy can be emitted 
as a V — V pair. Under the right conditions, this PBF process can be the dominant cooling agent 



in the evolution of a neutron star [)32[]. Such efficiency is due to the fact that the pairing phase 
transition is second order in nature. During the cooling of the star, the phase transition starts when 
the temperature T reaches when pairs begin to form, but thermal agitation will constantly induce 
the breaking of pairs with subsequent re-formation and possible neutrino pair emission. 




Figure 11: Left panel: Feynman diagram for v — V emission from the pair breaking and formation (PBF) 
process. Right panel: control functions Rpbf for the PBF process. 



^In the presence of a charged n condensate a new Urea neutrino emission pathway is open, see Table |l] and 



Eq (5.3). The development of a neutral condensate has, however, little effect on neutrino emission. 
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The emissivity of the PBF process can be written as 

^2 



™„ \lGim*rpf f 



T 1 ers -7 
3.51 X 10^' — I- X mfpFjTgafj Rj [Aj{T)/T] (5.6) 



cm^ s 



for a fermion / in a pairing state j ='S() or ^P2. The coefficients Ofj depend on the type of fermion 
and on vector and axial couplings Cy and Ca (see, e.g., [06|]). The control functions Rj are plotted 



in the right panel of Fig. |11[ These functions encapsulate the effect that the PBF process turns on 



when T reaches Tc and practically turns off at T ~ 0.2Tc when there is not enough thermal energy to 
break pairs. The PBF process was first discovered by Flowers, Ruderman, and Sutherland [ p3| ] and, 
independently, by Voskresensky and Senatorov [^1]]. It was, however, completely overlooked for 20 



years (?) until finally implemented in a cooling calculation in [44] and its importance emphasized 



in [32] 



An alternative way of looking at the PBF process is simply as an inter band transition of a 
nucleon [^]. Considering the particle spectrum in a paired state, right panel of Fig. ^ the lower 
branch (with e < £f — A) corresponds to paired particles whereas the upper branch to excited ones, 
i.e., coming from a "broken pair" which hence left a hole in the lower branch. A transition of 
a particle from the upper branch to a hole in the lower branch corresponds to the formation of a 
Cooper pair. 

6. Cooling of Superfluid Neutron Stars 



In order to complement the brief description of cooling presented in § |4.4{ , we now include 



the effects of pairing, but restrict ourself to cases in which no fast neutrino emission process is 



allowed (the minimal cooling scenario [|3^, ^). In the left panel of Fig. |12|, three examples of 
the thermal evolution of a 1.4 star are shown. The 4 cooling stages, introduced in Fig. ^, 
are similarly marked in this figure. The first model does not include any effect from pairing, and 
is very similar to the 1.4 Mq model of Fig. H (but not completely identical because of several 



improvements in the microphysics with respect to the older calculations of [33]). The other two 



models include nucleon pairing and in both cases the suppressing effects on the neutrino emission 



as well as the modification of the specific heat are taken into account (see Eqs. (5.3) and (5.2)), but 
the PBF process is artificially turned off in one case. When Cooper pair neutrino emission from 
the PBF process is turned off, one can clearly see the effect of the suppression of the modified 
Urea and bremsstrahlung processes during stage 3 which results in a much warmer star. One 
can also see the effect of the suppression of the specific heat which results in a faster cooling 
during the photon cooling era (stage 4). During stage 3, the specific heat is also suppressed, but its 
effect is not as pronounced as in stage 4. In the presence of the PBF process, however, the star is 
significantly colder during stage 3, the PBF process being more efficient than the MU process in 
the unsuppressed case (i.e., no pairing). In stage 4, the two trajectories with pairing, but with and 
without the PBF process, converge as cooling is now driven by the photons and only the specific 
heat suppression by pairing is relevant. Notice that during stage 2, the crust relaxation period, 
the model with pairing (and PBF) has a shorter relaxation time because of the suppression of the 
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Figure 12: Left panel: three models of cooling neutron stars illustrating the impact of pairing. Right panel: 
exploration of the effect of a late onset of neutron superfluidity The curves are labeled according 

to the maximum value of the neutron Tc'. Tc = max Tc{kF). No proton superconductivity is taken into 
account in these models. The star marks the observed value of for the Cas A neutron star at an age of 
330 yrs. See text for description. 



neutron specific heat from the ' So gap in the inner crust (the same effect is present in the model 
with pairing and no PBF, but is barely visible because the core is too hot). In stage L the tliree 
models give identical T/s as at that age Tg is controlled by the evolution of the upper layers where 
no superfluid is present. 

The constituent whose pairing has the major effect is the neutron fluid in the core. For il- 
lustrative purposes, the left panel of Fig. 12 employed a neutron gap chosen to maximize the 
differences between the three models: the bell shaped Tc{kf) curve (see the examples in Fig. 10) 
reaches a maximum of Tc = T^'^^ ~ 5 x 10^ K. The right panel of Fig. 12 shows a series of models 
in which the neutron gap Tc{kf) is rescaled by a constant factor, keeping the same shape. The 
dotted curves, with low values of Tc, exhibit a transit of the neutron star from the warm trajectory 
with Tc = toward the cold trajectories, with Tc = 1 or 2 x 10^ K. At early times when the star's 
core temperature T is above Tc, neutrons are normal but when T reaches Tc the phase transition 
starts at some location in the star (see the left panel of Fig. |T3| for a schematic of the evolution of 
Ly). At that stage, neutrons in a thick shell go through the phase transition and the neutrino lumi- 
nosity suddenly increases due to the triggering of the PBF process. The star then begins its transit 
toward a colder trajectory. As T decreases, this shell splits into two shells which slowly drift toward 
the lower and higher density regions away from the maximum of the bell-shaped Tc{kp) profile. 
When T is much below Tc{kp) everywhere in the core, Ly decreases rapidly. In the cases for which 
10^ K the onset of the phase transition occurs after the period of thermal relaxation of the 
crust and is seen as a second phase of rapid decrease of Te. For larger values of Tc, this transitory 
rapid cooling of the core is hidden by the crust which is not yet thermally relaxed with the core. 

The simple analytical solution of Eq. ( p^ ) gives some insight into this transitory behavior. 
When T > Tc, but <C To, the star follows the asymptotic "MU trajectory", T ~ t^^^^, and when 
T reaches Tc, at time t = tc, the neutrino luminosity suddenly increases. Despite the complicated 
T dependence of £™^, Eq. (5.6), the resulting luminosity, once integrated over the entire core 
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Figure 13: See text for description. Figure from [ p^ and 

(also, aided the bell shape of the Tc{kf) curve), is well approximated by a power law in the T 
regime in which some thick shell of neutrons is going through the phase transition. If we write 
L™P = / • L^u ^ fN^T^^ with / ~ 10, the solution of Eq. Q, replacing to by tc and To by Tc, 
gives 

^= [rq7(7^|:j7^ r = 109k(tmu//0'/' (when r«rc). (6.1) 

The central panel of Fig. [I^ shows the two asymptotic MU and PBF trajectories and the transit 
trajectory. Notice that in case the neutrino emission was already suppressed previously to the onset 
of neutron superfluidity, e.g., by proton superconductivity with a larger critical temperature, the 
MU trajectory is replaced by a warmer suppressed MU (SMU) trajectory and the initial part of the 



transit is much faster, as illustrated in the right panel of Fig. 12 



7. Cassiopeia A and its Cooling Neutron Star 

The Cassiopeia A supernova remnant (SNR) was discovered in radio in 1947 and is the second 
brightest radio source in the sky (after the Sun). It has since then been observed at almost all wave- 
lengths. Very likely, this supernova was observed by the first Astronomer Royal, John Flamsteed 
who, on August 16, 1680, when describing the stars in Cassiopeia constellation listed the star "3 
Cassiopeia" at a position almost coincident with the supernova remnant. This star had never been 
reported previously, and was never seen again until August 1999 when the first light observation of 
Chandra found a point source in the very center of the remnant (see Fig. [T^. 

The distance to the SNR is 3.4^q j kpc [p3|], and the direct observation, by the Hubble 
Space Telescope, of the remnant expansion implies a birth in the second half of the 17"' 



century [ ]14| ] and supports Flamsteed's observation. These observations give a present age of 331 
yrs for the neutron star in Cassiopeia A. The optical spectrum of the supernova has been observed 
through its light echo from scattering of the original light by a cloud of interstellar dust and shows 



the supernova was of type lib |24]. The progenitor was thus a red supergiant that had lost most of 
its hydrogen envelope, with an estimated zero age main sequence (ZAMS) mass of 16 to 20 Mq 
[ p2| , 11, ^] or even up to 25 Mq in the case of a binary system [56|. This implies a relatively 

1-11 > , A 1,^ HFTTi 



massive neutron star, i.e. likely ~ 1.4 [56|. 
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Figure 14: The Cassiopeia supernova remnant in X-rays: first light of Chandra, August 1999. (The 
neutron star is highhghted by the authors.) Image from ©NASA/CXC/SAO. 



The soft X-ray spectrum of the point source in the center of the SNR in Cassiopeia A is 



thermal, but its interpretation has been challenging [p9||. With a known distance, a measurement of 
the temperature implies a measurement of the star's radius, but spectral fits with a blackbody or a 
H atmosphere model resulted in an estimated radius of 0.5 and 2 km, respectively. It was only in 
2009 that a successful model was found: a non-magnetized^ C atmosphere, which implies a stellar 



radius between 8 to 18 km [22]. With this model, and analyzing 5 observations of the SNR, Heinke 
& Ho [ pO| ] found that the neutron star temperature had dropped by 4% from 2000 to 2009, from 
2.12 to 2.04 X 10^ K, and the observed flux had decreased by 21%. The neutron star in Cassiopeia 
A is thus the youngest known neutron star and the first one whose cooling has been observed in 
real time! 



8. Superfluid Neutrons in the Core of the Cassiopeia A Neutron Star 



The Chandra observations of Cassiopeia A presented in the previous section, together with 
its known distance, imply that the photon luminosity of the neutron star is 



lO^'^ergs-i . 



(8.1) 



With a measured ~ 2 x 10^ K we deduce an internal T ~ 4 x 10^ K from Eq. The 
star's total specific heat is thus C,, ~ 4 x 10^^ erg K^^ (from Fig. || or Eq. (O)). The observed 



ATg/Te ~ 4% [ pO| ] gives for the internal temperature AT/T ~ 8% over a ten years period since 
T ^ Tg. Assuming the observed cooling corresponds to a global cooling of the whole neutron star, 
its thermal energy loss is 



£th = Cyf ~ (4 X 10^^ erg K"') x (0.1 K s"') ~ 4 x 10" erg s 



37 



(8.2) 



There is, to date, no evidence for the presence of a magnetic field in the Cassiopeia A neutron star. 
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Figure 15: A very good fit of the rapid cooling of the neutron star in Cassiopeia A obtained assuming a 
recent onset of neutron superfluidity and the resulting increase in neutrino emission from the formation 
and breaking of Cooper pairs. The successful model assumes a maximum critical temperature 7c = 5.5 x 10^ 



K and the inset compares it with the six observational points, with their la error bars, from [201 and 
The two dotted curves with no neutron superfluidity, 7c = 0, and superfluidity with a higher 7c = 1 x 10^ K 
illustrate the sensitivity to Tc . Figure adapted from p7[ ] . 



which is 3-4 orders of magnitude larger than what it seen in Ly ! For a young neutron star, neutrinos 
are the prime candidates to induce such a large energy loss. 

The cooling rate of this neutron star is so large that it must be a transitory event, which was 
initiated only recently. Something critical occurred recently within this star! "Something critical" 
for a cooling neutron star points toward a critical temperature, and a phase transition is a good 
candidate. The results of the previous section exhibited a phase of accelerated cooling when the 
neutron pairing phase transition is triggered. With a Tc ~ 5 x 10^ K, a transitory cooling can 
occur at an age ~ 300 yrs as shown in the right panel of Fig. 12. 

This interpretation of the observed rapid cooling of the neutron star in Cassiopeia A as trig- 
gered by the recent onset of the neutron superfluid phase transition and the resulting increase in 
neutrino emission from the formation and breaking of pairs in the neutron superfluid was recently 
proposed in [37] and, independently, in |46]. Models such as the ones in the right panel of Fig. |T^ 
do not, however, exhibit a cooling rate as large as the observed one. 

A second key ingredient for reproducing the observed cooling rate is illustrated in the right 
panel of Fig. 13. The neutron star was very hot before the onset of neutron superfluidity. This is 
possible in the case the protons were already in a superconducting state, which implies that the 
corresponding critical temperature is significantly larger than the Tc for neutrons. That the for 
'So proton pairing is larger than the for neutron pairing is expected from the theoretical 
results presented in § 5J . However, a hot young neutron star can only be achieved in the case that 
the protons are superconducting in the entire core, so that neutrino emission form the modified Urea 
processes is strongly suppressed [46, 37]. This requirement places strong constraints on the proton 
'So pairing and is easier to fulfill if the neutron star mass is not too large as theoretical models 
show that proton superconductivity does not extend to very high densities. A better understanding 
of the progenitor of Cassiopeia A, and constraining the expected neutron star mass, is essential for 
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this scenario to work. 



'P2 



A very good fit to the observations is shown in Fig. 15, in which a maximum neutron 
pairing Tc of 5.5 x 10*^ K is employed and protons are assumed to be superconducting in the entire 
core with Tc{kp) > 10^ K everywhere. Very similar results were obtained independently in [46]. 



9. Discussion and Conclusions 

The observation, in real time, of the cooling of a neutron star is unique and its interpretation 
potentially imposes very strong constraints on the physics of ultra-dense matter. We have here 
given a basic presentation of the physical principles that are involved in trying to understand the 
interior of a neutron star. Neutrino emission processes and the likely occurrence of pairing are the 
two most important ingredients, and in both cases many unsettled issues still remain. 

We have presented arguments which, we hope, make it plausible that the observed 4% tem- 
perature drop, in a period of 10 years, of the young neutron star in the Cassiopeia A supernova 
remnant may signal the recent triggering of the phase transition of the core neutron fluid to a su- 
perfluid state. The resulting on-going formation of Cooper pairs results in a strong neutrino 
emission which can explain the observed rapid temperature drop of this neutron star. An essential 
requirement for the cooling to be as rapid as observed is that the neutron star was relatively hot 
at earlier times, and we described how this may be due to the suppression of the early neutrino 
emission by proton superconductivity. This would imply that the critical temperature, Tc, for pro- 
ton superconductivity is, everywhere in the star's core, larger than 10^ K. Under this condition, 
we found that the critical temperature for neutron superfluidity, which is density dependent, must 
have a maximum value of T^? ~ 5 x 10^ K. This would be the first direct evidence that neutron and 
proton superfluidity/superconductivity occur at supra-nuclear densities in the core of a neutron star 



[ ]37| , |46| ] and, further, these results would highly constrain their respective critical temperatures. 

Alternatively, one could consider the observed cooling of the star in Cassiopeia A to be instead 
due to a longer thermal relaxation timescale in some layer of the star than we have assumed. In 



such a case the estimate of Eq. ( |8.2D , which assumes an isothermal star, does not apply. As seen, 
e.g., in Fig. 12, the crust thermal relaxation occurs on a time-scale of a few decades. However, 
in case the crust thermal conductivity is smaller that the one employed here it may possible that 
this early temperature drop be delayed till the star is 300-400 years old (see, e.g., [^]). Similarly, 
the core thermal relaxation time may be much larger than usually considered: in case the inner 
core of the star cools rapidly it may take a few hundreds years for the star to become isothermal 
and this can induce a rapid decrease of Tg, as proposed in and presented in these proceedings 
by D. Blaschke. This scenario requires that the core thermal conductivity be lower than usually 
considered, by a factor of a few, and also implies that neutrons do not form a superfluid until the 
star is much colder. 

The neutron superfluid scenario we have presented fits within the minimal cooling scenario 



[35] and is compatible with observations of other cooling neutron stars [36]. The very similar 
model of [ p6| ] also successfully passes the same test as well as the opposite core relaxation scenario 
of [^. More work is required to confront these possibilities with other facets of the neutron star 
phenomenology. 
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